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Abstract: The Sakai-Sugimoto model provides a holographic description for chiral sym- 
metry breaking. We use this model to investigate chirally broken phases in an external 
magnetic field at finite isospin and baryon chemical potentials. The equations of motion 
for the bulk gauge fields are solved analytically and the free energy is computed from the 
Yang-Mills and Chern-Simons contributions to the D8 brane action. In the case of a neutral 
pion condensate, a magnetic field is found to induce nonzero gradients of the Goldstone 
boson fields corresponding to meson supercurrents. A charged pion condensate, on the 
other hand, expels the magnetic field due to the Meissner effect. Upon comparing the 
Gibbs free energies of these two phases we find that the rotation of the chiral condensate 
into a charged pion condensate for finite isospin chemical potentials is partially undone by 
switching on a magnetic field, and we determine the critical magnetic field which removes 
the charged pion condensate in a first-order phase transition. 
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1. Introduction 

Matter in large magnetic fields occurs in nature in the interior of neutron stars j^] , or can 
be created in the laboratory from noncentral collisions of heavy ions [0. In both contexts 
it is important to understand the response of strongly-interacting quark matter, ultimately 
described by Quantum Chromodynamics (QCD), to an external magnetic field. And in 
both contexts the interplay of the magnetic field with chiral properties of the matter, in 
particular with chiral symmetry breaking, are crucial. While heavy-ion collisions probe the 
regime of the chiral phase transition at large temperatures and small chemical potentials. 
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highly magnetized neutron stars ( "magnetars" ) are composed of dense and comparably cold 
matter, which also may be in a phase with broken chiral symmetry. In dense quark matter, 
chiral symmetry can be broken by the usual chiral condensates {ipip) or, in a three-flavor 
system at sufficiently large densities, by diquark condensates {tpip) in the color-flavor locked 
state Q. In this paper we shall consider a two-flavor system at flnite baryon and isospin 
chemical potential in the strong-coupling regime at large Nc, which may form different 
kinds of chiral condensates (i/^ip) depending on the values of temperature, the chemical 
potentials, and the external magnetic field. 

It has been shown using a Nambu-Jona-Lasinio (NJL) model that a magnetic field can 
act as a catalyst for chiral symmetry breaking |Q, |^, see also 0- Also chiral perturbation 
theory has been used to study the effect of magnetic fields ^, recently for instance in 
the context of the deconfinement and chiral phase transitions in Refs. [10| and |11, 12|, 



respectively. All these studies are restricted to the vacuum, i.e., they are done for the case 
of vanishing chemical potentials. Dense matter with nonvanishing chemical potentials in 
a magnetic field has been studied in the context of color superconductivity [l5|| , 
which, due to Goldstone boson currents and the axial anomaly, can be ferromagnetic |l6||. 



In this paper, we use the holographic model by Sakai and Sugimoto [17, at nonzero 
isospin and baryon chemical potentials to study the effect of a magnetic field on chirally 
broken phases.^ 

Holographic models have been used extensively to study the strong-coupling limit 



of gauge theories since the conjecture of the AdS/CFT correspondence ||2^, |2J], i.e., 
the duality between (the supergravity approximation of) string theories and (the strong- 
coupling limit of) conformal supersymmetric gauge theories; for a review see Ref. |2^]. The 
original and most prominent theory under investigation has been M = 4 super Yang-Mills 
theory which lives on the 3-1-1 dimensional boundary of AdSsxS^, and which is dual to 
type-IIB string theory living in this ten-dimensional space. While this supersymmetric 
gauge theory shares several properties with QCD, the differences to QCD are significant 
and thus it mainly serves as a model whose results should be compared with, not interpreted 
as, results from QCD. Therefore, a lot of effort has been spent to develop a gravity dual of 
QCD. Such a dual is unknown at present. Besides a "bottom-up" approach to AdS/QCD 



|2q , 27|, a promising model for capturing essential features of QCD is the Sakai-Sugimoto 



model, developed to study the physics of chiral symmetry breaking and meson spectra 



1 17, 18|. In contrast to the original AdS/CFT correspondence, but building upon an early 
proposal by Witten p8|, the Sakai-Sugimoto model deals with type-IIA string theory and 
a five-dimensional dual gauge theory, where by means of a compactified extra dimension 
supersymmetry is completely broken. This extra dimension separates massless fermions of 
right- and left-handed chirality, which are located on the intersections of D8 and D8 branes 
with Nc D4 branes. As all other known gravity duals, the model is dual to the large-A^c 
limit of the field theory, and the simple supergravity approximation, used in this paper as 
well as in most previous studies, corresponds to large values of the 't Hooft coupling. 



^For effects of magnetic fields in otlier holographic models of strongly coupled gauge theories with flavor 
degrees of freedom see e.g. Refs. [ p^ |2o| |2l| ]. 
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With nonzero chemical potentials and a magnetic field, we shall find meson super- 
currents and the Meissner effect in the chirally broken phases. Both phenomena are best 
understood as an analogy to (weak-coupling) superfluidity or superconductivity. For in- 
stance, a charged pion condensate of the form {d'y^u) can be viewed as Cooper pairing 
of two different fermion species, here an anti-down-quark and an up-quark. In general. 
Cooper pairing of two fermion species with chemical potentials /^i and ^2 takes place at 
a common Fermi surface given hy p, = (/Ui + ^2)/2. A mismatch in chemical potentials 
(5/i = {(ii — H2) /2 induces a "stress" on the pairing in trying to move the two Fermi surfaces 
apart. For not too large values of Sfi, the system can sustain the stress and the densities of 
the two fermion species are (at zero temperature) "locked" together, i.e., the difference in 
densities 6n = ni — 722 vanishes. For larger values of 5fi, and before completely breaking the 
condensate, the system may respond to the stress by leaving some, but not all, fermions 
around the Fermi surfaces unpaired, allowing for a nonzero 6n. The resulting state breaks 
rotational invariance, and it may even break translational invariance by giving rise to a 
crystalline structure. Anisotropic pion condensates in nuclear matter have been discussed 
a long time ago |29, |3^, crystalline structures of the superfluid order parameter are 

well-known in condensed matter physics |3^, 34] as well as in dense quark matter |35|, and 
also have been discussed in the context of chiral condensates |36], see also [^]. In either 
case, be it in a homogeneous manner or in a complicated crystalline structure, this un- 
conventional pairing induces nonzero "supercurrents" in the system, see for instance Refs. 
|38, 39, 40|. These supercurrents are cancelled by counter-propagating currents, typically 
coming from unpaired fermions, such that the net current in the system vanishes. 

In the case of a pion condensate of the form (dj^u), fl and 5fi correspond to the 
isospin, /i/, and baryon, hb, chemical potential, respectively. Consequently, one might 
expect anisotropic pairing upon increasing the "mismatch" fis- And, corresponding to the 
above 6n, a nonzero baryon number is expected. In the Sakai-Sugimoto model, a finite 
baryon number is taken into account via the Chern-Simons term. Localized baryons can be 
described by instantons of the effective gauge theory of the flavor branes pT| , 42, 43, 44, ^ 
corresponding to chiral skyrmions, which in the ground state form crystals rather than a 



liquid [14 6| , |47l , 48 1. However, we are interested in a homogeneous distribution of baryon 
(and isospin) density. It turns out that this can be achieved by a nonzero magnetic field 
in the model ||49| , 50 1, which is anyway of interest in the context of neutron star physics. 

A magnetic field, however, is expelled from the charged pion condensate because a 
condensate of charged bosons (be it Cooper pairs or, in our case, Goldstone bosons) acts as a 
superconductor and thus exhibits a Meissner effect.^ Accordingly, we shall find the above 
expectations of a supercurrent and nonzero baryon number not realized in the charged 
pion condensate which remains unmodified for (not too large) magnetic fields. A meson 
supercurrent as well as nonzero baryon (and isospin) numbers occur, for nonzero magnetic 
field, instead in the phase with a neutral pion condensate. This phase is energetically 
preferred over the charged pion condensate beyond a critical magnetic field which we shall 
compute. 

^Holographic models of superconductors and superfluids have recently been investigated in Refs. js^, 
[5^ , [5^ , [5^ , [ss] , see Ref. for a discussion of the Meissner effect. 
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Our paper is organized as follows. In Sec. § we briefly recapitulate the Sakai-Sugimoto 
model. We discuss our ansatz for solutions in the presence of baryon and isospin chemical 
potentials and a magnetic field and derive the equations of motion and the free energy for 
the chirally broken phases in Sec. |3[ The main part of the paper is Sec. ^. In this part we 



first discuss how to incorporate different chiral condensates into the model, see Sec. 4.1 



In Sec. 4.2 we solve the equations of motion for the sigma and the charged pion phase and 



compute their free energies. In particular, we discuss the Meissner effect in Sec. 4.2.2| . The 



results are used to discuss the currents and number densities in these phases in Sec. [4.3 
Finally, we compare their free energies to discuss the resulting phase diagram in Sec. 4.4 , 
and we give our conclusions in Sec. 

2. The model 

In this section, the basic equations of the Sakai-Sugimoto model shall be summarized. They 
will be needed in the subsequent sections. For more details about the setup of the model 
see for instance the original papers by Sakai and Sugimoto [17, p^j . 



2.1 Geometry of confined and deconfined phases 

The bulk background geometry is given by the ten-dimensional supergravity description of 
Nc D4 branes in type-IIA superstring theory compactified on a circle. There are two differ- 
ent solutions for the metric, realized in two different temperature regimes. The transition 
from one to the other is interpreted as the deconfinement phase transition. Similar to the 
original AdS/CFT setting at finite temperature |^], the deconfined phase has a black hole 
which is absent in the confined phase. 

The (euclidean) metric of the confined phase is given by p3] 



dslonf = (^) [dt"^ + Sijdx'dx^ + f{u)dxl] + ( - ) 



du 2 ,^2 

f{u) 



(2.1) 



Here, dVL\ is the metric of a four-sphere, and R is the curvature radius of the background 
which is related to the string coupling Qs and the string length via 

i?3 = TTQsNdl . (2.2) 

A crucial feature of the model is the compactified dimension X4 which has a radius which 
we can parametrize through the Kaluza-Klein mass 

3 u^"^ 

Mkk = - , Xi = Xi + 27r/MKK . (2.3) 

This breaks supersymmetry completely by giving Kaluza-Klein masses to the adjoint 
fermions of the dual gauge theory and the analogue of thermal masses to the adjoint 
scalars, leaving only gauge bosons in the spectrum of the low-energy limit as the latter 
are protected by gauge symmetry [28|. The point u = ukk is the tip of the cigar-shaped 
subspace spanned by X4 and the holographic coordinate n, and 

f{u) = 1 - ^ • (2.4) 
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The subspace spanned by the euchdean time t and the coordinate u is cyhnder-shaped, 
with the circumference given by the inverse temperature, t = t + 1/T. In the deconfined 
phase the coordinates t and interchange their roles, i.e., now the subspace spanned by 
X4 and u is cyhnder-shaped while the subspace spanned by t and u is cigar-shaped. In this 
case, the metric is 



u 3/2 / Ji\ 3/2 



RJ ' ' ' ' \uj 



+ u dO,i 



(2.5) 



where temperature is related to the tip of the cigar-shaped t-u space ut via 



4vr R3/2 

and 



(2.6) 



3 



f{u) ^ 1 - ^ . (2.7) 



The deconfinement phase transition is located at a critical temperature T = Tc where the 
free energies corresponding to the two phases are identical. This occurs at ukk = ut^ ^iid 
thus Tc = MKK/(2'7r). This critical temperature is independent of the chemical potential. 
Consequently, the model predicts a horizontal phase transition line in the T-fis plane, in 
accordance with expectations from QCD at infinite number of colors Nc [59]. 



The super gravity prescription depends on having the background weakly curved com- 
pared to the string scale. This is the case for large four-dimensional 't Hooft coupling 



A = g^^Nc = > 1 , (2.8) 



where the five-dimensional gauge coupling is given by g'^ = (27:)'^ ggig. 

The Kaluza-Klein mass sets the energy scale below which the dual field theory is ef- 
fectively four-dimensional. For large 't Hooft coupling, this scale is of the same order as 
the mass gap of the field theory; only for small A, where string corrections become impor- 
tant, does one have duality with non-supersymmetric large-A'^c QCD in four dimensions. 
However, there is already ample evidence that the limit of large 't Hooft coupling, where 
supergravity calculations are meaningful, does provide a useful tool for unravelling certain 
nonperturbative features of QCD. 

Sakai and Sugimoto |]l^ added Nj pairs of D8 and D8 branes which are transverse to 
the circle along X4. The intersections of these branes with the D4 branes carry massless 
fermions in the fundamental representation of the color group of opposite chirality, which 
are interpreted as massless quarks of QCD. As long as Nf <^ Nc, the D8/D8 branes can be 
treated as probe branes, i.e., the backreaction on the background geometry is neglected. 
Below we shall take Nf = 2. The D8 and D8 branes extend in all dimensions except for 
the coordinate X4 (whereas the D4 branes extend in the t, Xi, i = 1,...,4 directions). The 
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induced metrics on the probe branes in the confined and deconfined backgrounds are 



"•^DS.conf 



"*D8,deconf 



3/2 



du + u d^l^ 



r) 



f{u)dr + dijdx'dx^ 



3/2 



fin) 



■du + u d^4 



(2.9a) 



, (2.9b) 



where we abbreviated 



V{U) ^ ^1 + P{U) (|)'aX4)2 , V{U)^ yi+(|)'(5„X4) 



(2.10) 



Here the function X4^{u) gives the embedding of the D8 branes in the X/^-u subspace. 

The D4/D8-D8 setup provides the tools to study not only the deconfinement phase 
transition but also the chiral phase transition. In the X4 direction, the D8 branes are 
separated from the D8 branes by a distance L. The maximal separation of the branes is 
L = ir/MY^Yi in which case the branes are attached at opposite sides of the circle spanned 
by 2:4. Gauge fields on the D8 and D8 branes transforming under a local symmetry group 
U (Nf) induce a global symmetry group U{Nf) on the five-dimensional boundary at n = 00. 
More precisely, a gauge symmetry on the D8 branes induces a global symmetry at the four- 
dimensional subspace of the holographic boundary at X4 = 0, while the gauge symmetry 
on the D8 branes induces a separate global symmetry on the four-dimensional subspace 
at X4 = L. Therefore the total global symmetry can be interpreted as the chiral group 
U{Nf)L X U{Nf)R. 

So far we have viewed the gauge symmetry on the D8 branes as independent from that 
on the D8 branes. This is correct if the branes are geometrically separate. For example in 
the deconfined background, where the xa-u subspace is cylinder-shaped, the branes follow 
straight lines from n = up to u = 00, and thus are disconnected. However, it may also 
be energetically favored for the branes to be connected. In this case, the gauge symmetry 
reduces to joint rotations, given by the vectorial subgroup U{Nf)Lj^R. This is exactly the 
symmetry breaking pattern induced by a chiral condensate (for a discussion of the chiral 
condensate itself as an order parameter for chiral symmetry breaking within the Sakai- 



Sugimoto model see Refs. [ 6C , 31 1 ) . In fact, in the confined phase, where the x/^-u subspace is 
cigar-shaped, the branes must connect. In other words, chiral symmetry is always broken in 
the confined phase. Whether the branes are disconnected in the deconfined phase depends 
on the separation scale L. For sufficiently large L they are always disconnected, while 



for smaller L the connected phase may be favored for certain temperatures [32|. In other 
words, in the former case, deconfinement and the chiral phase transition are identical while 
in the latter case they differ and there exists a deconfined but chirally broken phase in the 
T-fj,B plane [^]. In this paper, we shall use maximally separated branes, i.e., L = tt/Mkk- 
This simplifies the treatment since in this case we always have duX4^ = because the D8 
branes follow geodesies. The case of not maximally separated planes, more precisely the 
limit where the radius of the compactified dimension is much larger than the separation 
distance, 1/Mkk ^ L, corresponds to an NJL model on the field theory side |64]. 
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Temperature and chemical potentials enter the model in very different ways. As ex- 
plained above, temperature has a geometric effect on the background metric, in particular 
a black hole forms for sufficiently large T. Chemical potentials, however, enter as boundary 
conditions for the gauge fields on the D8 and D8 branes, i.e., in the subsequent sections 
we will fix the baryon and isospin components of the temporal components of the "right- 
handed" and "left-handed" gauge fields at the boundary u = cxd by the isospin and baryon 
chemical potentials. Analogously, nontrivial boundary values of the spatial components 
of the gauge fields have the interpretation of spatial gradients in chiral condensates, cor- 
responding to supercurrents. We shall discuss the gauge field action associated with the 
flavor branes in more detail now. 



2.2 Yang-Mills and Chern-Simons action 

The total action for the D8 and D8 branes is given by the sum of the Dirac-Born-Infeld 
(DBI) and the Chern-Simons actions. As indicated in the introduction and as will become 
clear below, the Chern-Simons term is necessary to account for nonzero baryon and isospin 
numbers and, in our ansatz, will be proportional to the external magnetic field. For 
simplicity, we shall expand the DBI action for small gauge fields such that we obtain a 
Yang- Mills contribution instead. This was also done for instance in Ref. |5C], while other 



works used the full DBI action in a similar context, however for the simpler cases of a one- 



flavor system without isospin chemical potential |49| and without currents and magnetic 



field [65|. Our action takes the form 

'S'ds = 'S'ym + 5*03 • (2-11) 
Here, the Yang-Mills contribution is 

5ym = 2NfTsVi I (fx du e"* Vff (l " ^^^^^9^"" 9''"^AJ'vpJ'a^^^ , (2.12) 

where Tg = l/[(27r)^£3] is the D8-brane tension, where a' = l1, and where V4 = 87r^/3 is 
the volume of the unit 4-sphere. The remaining integrations are done over four-dimensional 
space-time t, xi, X2, 2:3, and over the holographic coordinate u. In the confined (deconfined) 
phase the limits for this integration are ukk{ut) < u < 00, and the factor 2 on the right- 
hand side of Eq. ( 2.12| ) accounts for integration over D8 and D8 branes. In this section. 



all expressions are thus valid for both confined and deconfined phase, which differ, besides 
the integration limits for u, by the metric g. The dilaton is e* = gs{u/R)^^'^, and the trace 
is taken over the internal U{2) space (from now on Nf = 2). Our convention for the field 
strength tensor is 

J' HP — d^j\i, di/j\^ , (2.13) 

where 1/ = 0, 1, 2, 3, u, and where is the U (2) gauge field. It is convenient to separate 
the U{1) part from the gauge fields and field strengths, 

A, = ^l + ^r^, :F,p = '-fl + ^T^, (2.14) 
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where a = 1,2,3 and are the Pauh matrices. With these conventions we have 
The Chern-Simons contribution in Eq. ( |2.11| ) is 



(2.15) 



'l27r2 
'967r2 



3 . 



1 



1 



ATr[F^] + -AF^ + -d 



iXr i2FA 



-A' 



d^x du \ ^i^ h^pF^y^ + \FupF„x 



+ 29, 



{ Ppa^t. + '^^abcApA'^A\ 



AiupcrX 



(2.16) 



where, in the first hne, we have used a notation in terms of differential forms in order to 
connect our expression to the one from Ref. |^] (our integration range is ukk < u < oo; 
therefore, in order to integrate over D8 and D8 branes we need an additional factor 2 in the 
prefactor compared to Eq. (2.8) in Ref. [^]). The change of numerical prefactors in going 
from the first to the second line comes from performing the trace and from our convention 
of the field strength ( 2.13| ) (the factors for the latter are hidden in the wedge products in 
the first line). 



3. Equations of motion and free energy in the chirally broken phase 

We can now derive the equations of motion for the gauge fields and the general form of 
the free energy, to be specified for various phases later. In this section and in the entire 
main part of the paper, we shall be concerned with the confined, i.e., chirally broken, phase 
whose metric g is given in Eq. (2.9a). For completeness we present the equations of motion 
and the free energy of the deconfined, i.e., chirally restored, phase in Appendix 

3.1 Equations of motion and ansatz including magnetic field, chemical poten- 
tials, and supercurrents 

We start by taking the variation with respect to the gauge fields of the Yang-Mills and 
Chern-Simons Lagrangians Cyu and Ccs ■ They are given by the integrands (including the 
prefactors outside the integral) of the actions in ( 2.12 ) and ( 2.16] ). We present the general 
form of the variations in Appendix Here we proceed by using those general expressions 
for our specific ansatz. 

The equations of motion obtained from the variations ( [A.2| ), ( |A.3| ), ( |A.6D , ( [A. 7 ) are 
complicated coupled nonlinear differential equations for the gauge fields. We shall now 
simplify these equations by transforming the holographic coordinate u, by choosing a par- 
ticular gauge, and by choosing a specific ansatz for the fields that captures the physics we 
are interested in. The new coordinate z we shall use from now on is defined through 



(uIk + UKKZ 



(3.1) 



We have z G [—00,00] while u G [ukk,oo]- In the new coordinate, the boundaries of the 
connected D8 and D8 branes correspond to z = —00 for X4 = ("left-handed fermions") 
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and z = +00 for x/^ = L = tt/Mkk ("right-handed fermions"), while the point z = 
corresponds to the tip of the cigar-shaped z-x/^ subspace in the bulk. We work in a gauge 



where Az = ^ |17, 1^]' U: ^^'^ ^ discussion of this choice. 

Now we specify our ansatz for the gauge fields. First, we set all components propor- 
tional to Ti and T2 in flavor space to zero and may then, for notational convenience drop 
the superscript 3 from the gauge fields and field strengths. Consequently, in the follow- 
ing we only have gauge fields and field strengths with a hat {A,F)^ corresponding to the 
1-components, and without any flavor index {A,F), corresponding to the Ta-components. 
This choice simplifies the calculations significantly but is a restriction for the possible chiral 



condensates we can capture, as we shall explain in Sec. [4.1 . 

The magnetic field is introduced as follows. The electromagnetic gauge group with 
generator Q = diag(gi,g2); where qi and q2 are the electric charges of the quark flavors, 
is a subgroup of U{2)l x U{2)b,- The magnetic field thus has baryon and isospin 
components, QBem = Bl + Bts, or 

B=^^B.^, B='^B.^. (3.2) 

We are interested in a system of up and down flavors, i.e., qi = 2/3 e, q2 = — l/3e with 
= 47r/137 and B = eBcm/Q, B = eBem/2, but mostly we shall derive general results, 
keeping B and B independent of each other. We should recall that the gauge symmetry 
in the bulk corresponds to a global symmetry at the boundary. Therefore, there is no 
electromagnetic gauge symmetry at the boundary, and in this sense Bem is not a dynamical 
magnetic field. 

We consider a spatially homogeneous magnetic field and, without loss of generality, let 
it point into the 3-direction. This requires nonzero field strengths F12 and F12. We can 
therefore choose the ansatz 

ii(x,z) = -xs^, i2(x,z)=xi^, (3.3a) 

^i(x,z) = -X2^, A2{^,z) =xi^, (3.3b) 

such that Fi2{z) = b{z), Fi2{z) = b{z), and the boundary values at z = ±00 of b{z), b{z) 
given by B, B. (Note that for non-constant b{z), b{z), we also have nonzero field strengths 

Fiz^ Fiz-) 

Next we account for the chemical potentials. This is done by relating the boundary 
values at z = ±00 for the gauge fields ^0(2;) and Ao{z) with the baryon and isospin chemical 
potentials /x b and fij p^, 33 , p3] . Consequently, we may have nonzero field strengths F^z , 



Foz- It turns out that within this ansatz nonzero values of the spatial gauge fields may 
be induced, i.e., we have to take into account A^lz), A3{z) and thus the field strengths 
Fsz, F^z- The boundary values at 2: = ±00 of the spatial gauge fields are identified with 
the gradients of the meson fields []T7| , |49t| . These gradients correspond, according to the 
usual hydrodynamic theory of a superfluid to "supercurrents" , i.e., currents of the 

condensate, in our context for instance the current of a pion condensate; see also Refs. 
1 38, Consequently, we shall identify (±00), ^3 (±00) with meson supercurrents j, 
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J. The supercurrents are not external parameters, hence we shall minimize the free energy 
with respect to them [4S, 50 1. They should not be confused with the "normal" currents 
Ji = (55efr/5^*, discussed in the Sakai-Sugimoto model in detail in Refs. 43 1, and 
computed below in Sec. 4.2.1. The supercurrents rather act as a source for the normal 
currents. 

Now we can insert the ansatz and the coordinate transformation ( |3.1| ) into the general 
equations of motion ( [A.2| ), ( |A.3| ), ( [A.(^ ), ( [A. 7 ). We then need to replace — > iAo, since 
we are working in euclidean space. We find the following equations for the magnetic field 



dMz)dzb] = dMz)dzb] = Q., 



where 



k{z) 



(3.4) 



(3.5) 



They arise from the Yang-Mills variation with respect to the spatial gauge field, Eqs. 
( |A.2b| ) and (A. 3b) and contain no contribution from the Chern-Simons term. Moreover, 
they decouple from the equations for the other fields, which are 



au 



KK 



au 



KK 



au 



KK 



au 



KK 



where 



a = 



b{z)F,3 + b{z)F,3 
b{z)F,3 + b{z)F,3 
b{z)F,o + b{z)F,Q 
b{z)F,o + b{z)F,o 

277r 
2A"' 



(3.6a) 
(3.6b) 
(3.6c) 
(3.6d) 

(3.7) 



In all four equations in ( |3.6D the left-hand side comes from the variation of the Yang-Mills 
contribution, while the right-hand side originates from the Chern-Simons contribution. We 
see that the latter is proportional to the magnetic field. The equations ( ^^ ) and ( |3.6D shall 
be solved analytically in Sec. Before doing so we use these equations to derive a simple 
expression for the free energy. 

3.2 Free energy and holographic renormalization 



With the metric of the confined phase ( |2.lD , the relations between the parameters of the 
model ( |2.2D , ( ^ ) , ( ^ ) , and the new coordinate z ( |3.1D , the Yang- Mills part of the action 
( 2.12D can be written as 



dz 



16M2j^A;2/3(z) 



+ 



9(27rQ')2nKK ^kk 
where = 0, 1, 2, 3. Here, k{z) is given in Eq. (p. 5]), 



'^'^k{z)TT[J^%] + -h{z)TT[T%] 



h{z) = (uIk + UKKZ'^) 



-1/3 



(3.^ 



(3.9) 
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and 



2167r3 



(3.10) 



In deriving Eq. ( |3.8D we have used that the field strengths are symmetric or antisymmetric 
functions of z. We shall see later that this is indeed the case for all phases we consider. 
This form of the action is general, and it is straightforward to insert our ansatz discussed 
in the previous subsection. 

To compute the Chern-Simons contribution to the free energy we first note that the 
surface term (last term on the right-hand side of Eq. ( p.l6|) ) gives a nonzero contribution. 
Within our ansatz the term oc (i(y4Tr[j43]) vanishes since our only nonzero flavor compo- 
nents of the gauge fields are proportional to 1 and t^; however, the term oc (i(74Tr[FA]) 
does not vanish. We find 



a^(i,Fp,AA)e^'^^'^^ = 2b{AsF,o - A^F,^ + 2AqF,^ - 2A^F,^) 

+ 2b{A3F,o-AoF,s). 



Inserting these expressions into the Chern-Simons action ( |2.16| ) yields, with 

iVc/(167r2) = OK, 



(3.11a) 
(3.11b) 

(3.11c) 
iAo and 



OK 



5*03 = / dx^ 



kM^V 
34k T 



dz 



b i^AsF^o + AsF^o - AoF,3 - AqF.^ 
+ b (isF^o + A3FZ0 - AoFz3 - AoFz3 
dzk{z)iF^, + F^,-F^3-F^3) 

k{z){AoF,o + AoF,o - A3F,3 - ^gF^g 



z=+oo 
z=—oo 



(3.12) 



where, in the second step, we have used the equations of motion (|3.q ), and where V is the 
three-dimensional volume of space and T the temperature. In changing the integration 
over the holographic coordinate from u E [ukkiOo] to z £ [— oo,cxd] we have assumed that 
the integrand is symmetric in z. In all phases we consider this turns out to be the case. 
Putting the Yang-Mills and Chern-Simons contribution together, we obtain the free energy 
density O = T{Sym + Scs)/V 



= +f^,+ ''^^^KK 



37/2 



dzk{z){-Fl - Fl + F^,, + Fl) 

z=+oo 



k{z){AoFzo + AoFzo - A3FZ3 - A3FZ3) 



(3.13) 



where the geometric contribution is given by the field-independent first term on the 
right-hand side of Eq. (|3.8D . This term is independent of all gauge fields and field strengths 
and thus plays no role in discussing the physical properties of a given phase. Moreover, we 
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shall only compare free energies of phases with identical embedding of the flavor branes. 
Hence, for our purpose, this term can simply be dropped from now on. The term il.h in 



Eq. ( 3.13 ) is given by 



+ 



dz h{z)[b\z) + b\z)] 



4u^j^ / dxidx2 



dzkiz)[{d,bf + id,bY 



(3.14) 



Both contributions of 0,^ solely depend on the magnetic field (remember that the equations 
of motion for b and b (3.4) decouple from the other field equations). Therefore, Qb is 
irrelevant for minimizing the free energy with respect to the supercurrents j and j. However, 
it can play a role when comparing free energies. This poses a problem, as both terms of 
Qb are divergent. 

Let us first consider only constant functions b{z) = B and b{z) = B, for which only 
the first term in Qb is present. Since we have already divided by the volume V of 3-space, 
we would expect a finite energy density from a homogeneous magnetic field, but because 
of the extra holographic dimension, this is not the case. In fact, since h(z) ~ the 
divergence of fi^ comes from the jz] — > 00 limits of integration and is thus a typical holo- 
graphic divergence which can be treated by holographic renormalization [^0[. Here we do 
not attempt to provide a complete discussion of this procedure, which for the (nonconfor- 



mal) Sakai-Sugimoto model has been introduced only recently |71, 72]. We rather follow 
the method outlined in these papers and subtract a counterterm, fixed by a physical renor- 
malization condition, as follows. After restricting the holographic integration in Qb to a 
finite interval —A < 2: < A, we subtract a counterterm 50;, (A) which cancels the divergence 
and obtain a renormalized contribution r^jj^"^. We also include a finite counterterm which 
is fixed by requiring the free energy in the absence of any chemical potential to vanish. 



0(/iB,/ = 0) =0. 



(3.15) 



This condition is motivated by the observation that 0, should be the matter part of the free 
energy, i.e., it should describe the fermions and their interaction with the magnetic field. 
In particular, we thus require that the energy density of the (nondynamical) magnetic field 
in the absence of any matter be left out. This we shall later treat separately when we 
consider the Gibbs free energy (the Legendre transform from fixed internal magnetic field 



to fixed external magnetic field) in Sec. iA. The condition ( 3.15 ) implies that we have to 
require 

= = lim [QbW - Snb{A)] . (3.16) 

A^oo 



To find the exact form of the counterterms we first note that, for constant b{z) 
b{z)=B, 

AV3 _r(5/6) 



QbW 



3k(^2 ^ ^2^ 



u 



1/3 
KK 



r(i/3) 



+ 



"KK 
A5/3 



B and 



(3.17) 



The counterterm (5 fit (A) should depend only on fields and geometric data on the slice 
z = A, in particular it should only involve the induced metric 7^jy on the slice and not 
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the complete metric g. By including appropriate factors of the dilaton 72] and an 
appropriate numerical factor to fulfill the condition ( |3.16| ) we find 



5Qb(A) 



R 
2" 



<I.\ 1/3 



7rr(5/6) /UKKy/2 

r(i/3) V R ) 



<I.\ -1/3- 



C(A) 



(3.18) 



with 



C(A) 



T^Vi{2'Ka') 



l\2 



(Jflj 



AV6 

1/6 



+ 



u 



KK 



All/6 



(3.19) 



where we have used Eq. ( A.4| ) and where the indices ix,v,p,a run over 0,1,2,3. With 
this counterterm, the term proportional to A^/^ (A*^) in Eq. ( p. 17 ) is cancelled by the first 
(second) term in Eq. ( |3.18| ) 

In the case of a magnetic field which is not constant in the bulk, the second term in 
Ob as given by Eq. ( 3.14D is also divergent, but its divergence comes from the integration 
over the spatial directions perpendicular to the magnetic field, regardless of whether the 
holographic z-integration is finite or not. Therefore, we cannot treat this term by the usual 
holographic renormalization and we interpret this divergence, when present, as a Meissner 
effect: a phase where a homogeneous magnetic field B^ra, which fixes the boundary values 
of h{z) and b{z)^ is only possible for non-constant functions in z, is infinitely penalized 
such that only Bem = is allowed. As we shall see, this will be the case for the charged 



pion condensate, to be discussed further in Sec. 4.2.2. At this point we already observe 



that the role of the spatial directions transverse to the magnetic field is no coincidence. It 
points to the necessity of currents in these directions which produce a magnetic field equal 
in magnitude but with opposite direction compared to the external magnetic field. This 
leads to a vanishing total magnetic field in the system, which is nothing but the Meissner 
effect for superconductors. 



4. Chirally broken phases in a magnetic field 

In this section we solve the equations of motion for the chirally broken phase. We shall 
distinguish between two different chirally broken phases, the a and the vr phase. This is 



the main part of the paper, and the main physical results can be found in Sees. 4.3 and 



4.4 



4.1 Chiral rotations and resulting boundary conditions 

In Nf = 2 chiral perturbation theory the chiral field U £ U{2) describing the Goldstone 
bosons is given by 



U 



JiV+'Para)/fTT 



(4.1) 



where f-j^ is the pion decay constant (in the Sakai-Sugimoto model, /^r = 2Mkk\/'^/'?i" 
1 17, |4^). The r] meson (the r/' for Nj = 3) becomes massive in QCD due to the explicit 
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breaking of the U{1)a through the axial anomaly. This is realized in the Sakai-Sugimoto 
model through the Chern-Simons term, and the mass of the rj can be computed within the 
model, = AMkk V^/ /Nc/{3V3tt) 0, see also Refs. 0, |7|, 



In the Sakai-Sugimoto model the chiral field is given by the holonomy |17] 



oo 



U = Pexp\^i J dzAzj . (4.2) 

As mentioned above, we work in a gauge where Az = 0. This is only possible by using the 
full U{2)l X U{2)fi symmetry, implicitly taking into account a 9 term, see also comment 
below Eq. ( |4.(]| ). It seems we can then only consider the vacuum U = 1. However, we can 
keep the Az = gauge and recover other vacua encoded in the boundary values of the 



gauge fields. This is explained in detail for instance in Ref. [66|. We shall now recapitulate 
this explanation and apply it to our case. 

Consider a potential V[ij,l, Hr, U ((/>)] which is invariant under U {Nf)^ x U {Nf)fi. Here, 
f^L, ^J'R G U{Nf) are fixed external parameters. For the following argument we denote these 
parameters simply by fiL, reminiscent of the chemical potentials, but one should keep 
in mind that this notation also includes the magnetic field. The chiral field U is written as 
a function of a parameter (p with respect to which we have to minimize the potential to find 
the vacuum. This parameter is a symbol for the meson fields in Eq. ( [4.1[ ). The external 
parameters transform under the global symmetry as fii — > gJ^^fJ-LdL, ^i? dR^fJ-RdR, while 
the chiral field transforms asU ^ dJ^UgR, where gi G U {Nf)L, gji E U {Nf)ji. Via a global 
symmetry transformation we have V[hl, Hr,U {(/))] = V[gl^{4))iJLLgL{4>)ig'^{4>)l^RgR{4>): 1] 
with 0-dependent transformations gL{4>), gR{4') such that g^^{(t>)U{(j))gR{(l)) = 1. To find 
the vacuum of the theory it obviously does not matter whether we use the original poten- 
tial or the potential with the transformed quantities because both expressions are simply 
identical. Consequently, instead of keeping the external parameters fixed and varying the 
chiral field we can fix the chiral field to be the unit matrix and vary the external param- 
eters. Of course we cannot simply treat the external parameters as arbitrary continuous 
quantities with respect to which we minimize the potential. We need to ensure that they 
are connected by a transformation to their physical values. We shall see below that within 
our ansatz the allowed rotated parameters only assume two discrete values, such that we 
simply have to compare two separate phases with each other. After minimization of the 
potential, the physical vacuum is given by applying the rotation found from minimization 
"backwards" onto the unit matrix, i.e., 

U = gLg-R^ ■ (4.3) 
Without loss of generality we can set gR = \ and thus U = gi- We can write 

= e'^V+'fiaTa)/h = (.iv/h ^ + ^^""^^ ^ (4,4) 

where u//^ = cos(</j//^), vTa//^ = yja/v? sin(99//^) with ip = {ipl + Lpl + ^plY^'^. This is the 
usual form of the chiral field in chiral perturbation theory, where the massive mode, the 
"sigma" , is frozen and the effective theory describes the remaining meson modes. Therefore, 
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io(±oo) 


Ao(±oo) 


b{±oo) 


b{±oo) 


i3(±00) 


^3(±oo) 


a 






B 


B 


±2i 




vr 




±2nj 


B{0) 


±B (0) 


±2j 






Table 1: Boundary conditions in the sigma and pion phases. The boundary conditions for the 
temporal components of the gauge fields correspond to the baryon and isospin chemical potentials, 
while the boundary conditions for the field strengths F12 = b, F12 = b correspond to the baryon and 
isospin components of the magnetic field. The boundary conditions for the spatial components A3, 
are given by the meson supercurrents j, j. These currents are not external parameters but have 
to be determined by minimizing the free energy. In the tt phase A3 (±00) has to vanish to ensure a 



well-defined behavior of the gauge fields under parity transformations, see Eq. (C.6) and discussion 
above this equation. For a discussion of the normalization of the chemical potentials see Eq. ( 4.34 ) 
and below. The zeros in parantheses for the magnetic fields in the charged pion condensed phase 
indicate that eventually we shall set B = B = because of the Meissner effect in this phase, see 



Sec. 4.2.2 



both sides of Eq. (4.4) contain four degrees of freedom; for the right-hand side we have the 
condition cr^ + vr^ = which is obvious from the definitions of a and TTa- 

To apply a rotation given by on the (left-handed) external parameter ^u^ note that 
our physical chemical potentials and the magnetic field are diagonal in flavor space and 
identical for L and R, fii = IJ-R = A'sl + IJ-iTa, Bem,L = Bem,R = Bl + Bt^. The baryon 
part oc 1 does obviously not change under a U{2) transformation. We thus only have to 
consider how the isospin part tx transforms. We find 



QL^T^igL = [7r+(7r° + ia)T+ + tt (vr" - ic7)r_ + (1 

Jtt 



27r+^-)r3 



(4.5) 



where r± = ri it iT2, and where we have introduced the neutral pion vr'' = tts and the 
charged pions vr^ = vri =p i7r2. In our ansatz described in Sec. 3.1 we have restricted 
ourselves to diagonal gauge fields. Since the chemical potentials and the magnetic field are 
the boundary values for the gauge fields, they have to be diagonal too. Consequently, we 
can only allow for transformations ( |4.5D that transform into a matrix oc . There are two 
(nontrivial) possibilities to make the coefficients in front of t_|_, r_ vanish: (i) vr"*" = tt" = 
which leads to gJ^T^QL = T3 and (u) vr^ = o" = which leads to g^^r^gL = —T3- Hence 
we can either leave the isospin components of the chemical potentials and the magnetic 
field invariant or flip their sign. This means that the parameter <j) in U{(j)) above is in fact 
discrete, not continuous. Had we allowed for off-diagonal components in the gauge fields, 
we could have described arbitrary linear combinations of the pion fields. 

These somewhat formal arguments have a very intuitive geometric interpretation |^5| : 
another (simpler, but less precise) way of saying what we have just explained is the fol- 
lowing. Think of the D8 branes as a left-handed up-brane and a left-handed down-brane 
and of the D8 branes as a right-handed up-brane and a right-handed down-brane. Then, 
a chirally broken phase can be constructed by connecting (i) the left-handed up-brane 
with the right-handed up-brane and likewise for the down-branes or (ii) the left-handed 
up-brane with the right-handed down-brane and vice versa. These two possibilities corre- 
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spond exactly to the two cases from the above formal argument: case (i) corresponds to a 
condensate where equal quark flavors participate, i.e., a combination of a and vr*^ with u — u 
and d — d pairing. In the remainder of the paper we shall refer to this case as the a phase. 
Case (ii) corresponds to a charged pion condensate with nonzero (d'j^u), {u^^d), to which 
we shall refer as the vr phase (this phase is sometimes called "p" |^, ) . As a summary of 
this section and a reminder for the subsequent sections, we present the resulting boundary 
conditions for the a and the tt phases in Table |^. 

In this Table we have also included the supercurrents j, j, which, in our gauge ql = U, 
have the form gJ^VgL 0- With Eq. (g^) this becomes for the two phases 



(i) a phase : -igl^VgL = ^ + ^ (f^V^° - 7r°Va) , (4.6a) 



(ii) vr phase : -i9L^^9L = + {tt Vvr+ - vr+Vvr ) 



(4.6b) 



We see that the supercurrents are diagonal, i.e., our ansatz with nonvanishing 1-component 
j and T3-component j is consistent. Interestingly, an anisotropic rj condensate appears in 
the 1-components j. The tj condensate has dropped out in Eq. (^?5|), and thus our boundary 
conditions, given by the chemical potentials and the magnetic field modified by the rotation 
( |4.5| ), do not reveal whether there is an admixture of an rj condensate in the a phase. On the 
other hand, a nonzero supercurrent j seems to indicate the presence of an rj supercurrent. 
Indeed, we shall see later that in the a phase a nonzero j is induced. The term Vr/ in 
Eqs. ( [4.6D appears due to our use of the full U{2)i x U{2)ji symmetry. Strictly speaking, 
our Lagrangian breaks the axial U{1)a because of the presence of the Chern-Simons term. 
However, this symmetry is preserved if one compensates a C/(l)^ rotation by a shift of the 6 
parameter, whose realization in the Sakai-Sugimoto model is discussed in Ref. |17, 76|, 



see also Ref. |78]. We thus implicitly adjust the 6 parameter when using the full gauge 
symmetry, absorbing a constant rj mode into 6. We shall proceed within this simplification, 
but have to keep in mind that in a more complete approach one would have to consider a 
fixed 9 and allow for a constant r] mode explicitly. Such an approach would be of interest 
especially in view of recent studies of possible (CP-violating) ry condensates in an NJL 



model calculation [79|, or, including a magnetic field, in the linear sigma model [12|. 

We finally remark that our setup does not include the possibility of diquark condensa- 
tion of the form (ud) , which is expected to lead to color superconductivity of quark matter 



at sufficiently large baryon chemical potential |8C]. However, color superconductivity does 
not necessarily occur in the large Nc limit where a "chiral density wave" is a strong candi- 
date for the ground state or, as suggested recently, quark matter may be confined 
even for large chemical potentials |5^, 83]. 



4.2 Solutions of the equations of motion and free energies 



We can now solve the equations of motion (|3.4| ) and (|3.6D for the two sets of boundary 
conditions given in Table ||. For notational convenience we set ukk = 1 (the final results 



in Sees. 4.3 and 4.4 do not depend on ukk and thus all physical quantities will have the 
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correct dimensions) . For both sets of boundary conditions we first note that the differential 
equations ( |3.6| ) can be solved by defining the new functions 



Ft{z)^k{z)^:^^^. (4.7) 



Then, the four equations (|3.6|) are equivalent to 



± _ a[b{z) ± b{z)] , n F^- "[K^) ± Kz)] j,±, . . 

Now the two equations with the upper sign are decoupled from the two equations with the 
lower sign. To proceed, we have to distinguish between the two chirally broken phases. 

4.2.1 Sigma phase 

With the boundary conditions of the a phase from Table |^ and with Eqs. (3.4) we conclude 



that the magnetic fields are constant in the bulk 

Kz) = B, b{z) = B. (4.9) 
In the following, we shall denote the dimensionless magnetic fields by 



We can now solve Eqs. ( |4.8| ) for completely general boundary conditions for the gauge 
fields. This is done in Appendix ^, where we present some technical details. Here we 
proceed with the specific solution obtained from the boundary conditions given in the first 
row of Table This solution yields the gauge fields 

Ao{z) = +i[C+(z) + C.{z) - T+]+j[C+{z) - C_(z) - r_] , (4.11a) 

^o(^) = 2^iI + J[C+{z) + C.{z) - T+]+j[C+{z) - C^{z) - r_] , (4.11b) 

i3(^) =3[S+{z) + S-{z)]+iS+{z)-S,{z)], (4.11c) 

A:,{z)=j[S+{z) + S-{z)]+3[S^{z)-S^{z)], (4.11d) 

where we have abbreviated 

_ cosh[(i? lb i?) arctan z] q ( \_ sinh[(i? it i?) arctan z] 

^±\Z) = ~ 1 '-'±{^) = ~ 1 (4.iza) 

sinh[7r(S±S)/2] sinh[7r(5 ± B)/2] 

, tt(B + B) tt(B-B) 
T± = coth ^ ^ ' ± coth — . (4.12b) 

Note that T± = C+(cxd) ib C_(oo) = C+(— oo) it C_(— oo). Since the functions C±{z) and 
S±{z) are symmetric and antisymmetric in z, respectively, both temporal components of 
the gauge fields are symmetric while both spatial components are antisymmetric. Together 
with the behavior of the supercurrents under a parity transformation this ensures that the 
gauge fields transform as a vector under parity, see discussion below Eq. ( [B.9[ ). We plot 
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Figure 1: Energetically preferred configuration of the gauge fields as a function of the holographic 
coordinate z for the sigma phase (left panel) and the charged pion phase (right panel). For the 
sigma phase we have chosen a dimensionless magnetic field eBcm = 2. In the pion phase, i?om = 
due to the Meissner effect. The boundary values for Ao{z) and Ao{z) are given by (twice) the 
baryon and isospin chemical potentials, respectively. The boundary values of ^3(2) and ^3(2;) yield 
the meson supercurrents and are determined dynamically from minimization of the free energy. 



the gauge fields with the supercurrents determined from minimization of the free energy, 
see Eqs. ( |4.15D , in the left panel of Fig. ||. 

Next, we insert the gauge fields and the resulting field strengths into the free energy 



(3.13). We drop the contributions Jig and, via holographic renormalization, Of,, as explained 
in Sec. |3.2| . Then we obtain (for details see Appendix ^ 



n 



4pB{jB + jB) - 4pi{jB + jB)] , 



with 



2smh^[-K{B±B)/2] 



(4.13) 



(4.14) 



The asymptotic values of the functions p±{B,B) at small and large magnetic fields are 
shown in Table ^ in Appendix ^ 

Minimizing SI with respect to j, j yields 



PB + Pi B + B Pb - P-i B-B 
3 = + 



2 p+{B,B) ' 2 p_{B,B)' 
PB + Pi B + B pb - Pi B-B 



2 p+{B,B) 



2 p^{B,B) 



(4.15a) 
(4.15b) 



One can check that this is indeed a minimum of Vt: the matrix of second derivatives of 
O with respect to the supercurrents has eigenvalues 8kM;^j^p+/3, 8kM^j^/3_/3, which are 
independent of j and j and positive for all B, B. 
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As already mentioned below Eq. (|3.3| ) we recall that the supercurrents j, j act as a 
source for the normal currents J , J which are the spatial 3-components of the four-currents 



where the Yang-Mills and Chern-Simons contributions are given by p^ , 

Nr. 



^YM,R/L 



=±oo 



7-M 

^CS, R/L 



247r2 



I/-' per 



=±oo 



(4.16) 

(4.17a) 
(4.17b) 



Here, the indices iJ:,u,p,a run over 0,1,2,3, the upper (lower) signs correspond to R (L), 
and we have, in the Chern-Simons contribution, already used that in our ansatz the off- 
diagonal components of the gauge fields in flavor space vanish. With the gauge fields ( 4.11| ) 
and the field strengths ( p.9D we obtain the baryon and isospin components of the spatial 
currents, 



J = Jr = -Jl = kM. 



KK 



' ■ coth 



P+{B,B 
B) 2 



P-{B,B) 



(4.18a) 



J = Jr 



-Jl 



[fiB + Pi) — coth ■ 



P4 



B,B) 



{PB - PI, 



{B 



B? , tt{B 
' ■ coth — 



B) 



-{ixbB + ixiB) 



(4.18b) 



P-{B,B) 

where the terms with prefactor 2/3 are the Chern-Simons contributions. These currents are 
already evaluated at the minimum of the free energy, i.e., we have inserted the supercurrents 
( |4.15 ). They add up to zero in the sums Jl + Jr and Jl + Jr, corresponding to vanishing 
baryon and isospin currents, however they yield nonzero axial currents. 

For small magnetic fields, the linear terms of the Yang-Mills and Chern-Simons con- 
tributions cancel exactly and thus the currents become cubic in the magnetic field, 

^2 



J 



J ^ 



54 
54 



{B + BYifiB + pi) + {B- BYifiB - pi) 
{B + B)^ifiB + pi) -{B- B)\fiB - pi) 



(4.19a) 
(4.19b) 



The cancellation of the linear terms seems to suggest that there might also be cancellations 
in the cubic terms from terms we have neglected upon expanding the DBI action for small 
gauge fields, possibly leading to vanishing currents. By considering a one-flavor system 
(where the use of the full DBI action is much simpler) we have checked that this is not the 



case." 



^This is at variance with Ref. ^] where a vanishing axial current has been found. In this reference, 
however, a modified action has been used where certain surface terms are dropped, essentially to force the 
current to vanish. Since we do not see any contradiction in a nonvanishing current (in the presence of an 
external magnetic field) we leave the full resolution of this discrepancy to future studies. 
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We plot the currents J ^ J in Fig. ^, where we also show the supercurrents and the 
densities in the sigma phase. We remark that the expansion of the 3-component of the 
gauge fields does not contain the complete current. One rather only finds the Yang-Mills 
contribution as a coefficient in front of the next-to-leading term of the asymptotic expan- 
sion, 

Uz) = ±2i^ %^ - + o(^), Mz) = ±2j^^^ - + o(^). (4.20) 



Finally, inserting the values ( [4.15 ) back into yields the value of the free energy at 
the minimum. 



^ 2kM|k 



(4.21) 



We see that for vanishing magnetic fields Vt^j = 0, i.e., the free energy does not depend on 
any of the chemical potentials. This is the expected result for the sigma phase and has 
also been observed in Ref. [ |6^ . 

4.2.2 Pion phase and Meissner effect 

In this case the boundary conditions are given in the second row of Table ffl, and the 



differential equations (3^) for the magnetic fields have the solution 



2B 

h{z) = B , b{z) = — arctan z . (4.22) 

TT 



As we have discussed at the end of Sec. 3.2, nonconstant functions 6 or 6 lead to an infinite 
contribution to the free energy which cannot be removed by holographic renormalization, 
but which enforces a vanishing magnetic field, indicating a Meissner effect. In the vr phase, it 
is the nonconstant isospin component h{z) which leads to this conclusion. This is only to be 
expected since the condensate of pions carries an electric charge, and thus the system is an 
electromagnetic superconductor. By the Meissner effect, a magnetic field is induced which is 
opposite, but equal in magnitude, to the applied magnetic field, such that B^m = and thus 
B = B = 0. Of course our electromagnetic group is only global and thus the microscopic 
description of the Meissner effect, for instance in terms of a Meissner mass for the photon, 
is not straightforward. However, in terms of supercurrents, the effect can be described 
quite naturally: in fact we have to allow for a supercurrent in the directions transverse 
to the magnetic field, i.e., j^(x, z) = ^6(z) (x2, — xi, 0), such that curlj^ = —b (and the 
same for the components jg, b). This is the usual London equation for a superconductor, 
see for instance Ref. |84]. Consequently, we need to add the supercurrents jg, to the 
boundary conditions of the gauge fields Ai(x, z), A2(x, -z), Ai(x, z), 742(x, z) from Eqs. 
( |3.3| ) such that the total boundary conditions (and thereby the total magnetic field in the 
superconductor) vanish, B = B = 0. This condition renders the equations of motion for 
the pion phase very simple. We shall, however, solve these equations for arbitrary magnetic 
fields and only at the end set B = B = 0. This provides us with a better understanding of 
the structure of the solution, for instance its behavior under parity transformations. 
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We defer all technical details and the solution for general boundary conditions with the 
magnetic fields ( [4.22| ) to Appendix]^. For the specific boundary conditions characterizing 
the charged pion condensate we find the solutions 

Ao{z) = 2f,B + MC+{z) + C-{z) - f+] +j[C+{z) - C-{z) - T_] , (4.23a) 

Ao{z) = fii[S+{z) + S^{z)] +j[S+{z) - S^{z)] , (4.23b) 

A^iz) =j[S+{z) + S^{z)]+fij[S+{z)-S^{z)], (4.23c) 

A^{z) = j[C+{z) + (7_(z) - f_] + fii[C+{z) - C-{z) - f+] , (4.23d) 

where we abbreviated 

5, (.) . . (.) . «±4^ei) , (4,24b) 

with 

^ TT TTsi / vr^ ± 25 arctan z \ /. „^ n 

"^^f^' " Vl' " i i73 j ■ 

F,(=) . f -^^^^"-' ) , (4,25b) 

and = (5+(ibcxD), = P+(iboo). Here, erf is the error function and erfi(z) = erf(iz)/i. 
The functions C± , S±, T± are the more complicated counterparts of the functions C± , S± , 
T± from Eqs. ( [4.12 ). They share the same property T± = Cj^{oo) ± C-{oo) = C'+(— oo) it 



C_(— oo), and, as their counterparts, C±{z) and S±{z) are symmetric and antisymmetric in 
z, respectively. This means that the temporal 1-component and the spatial rs-component 
are symmetric in z while the temporal rs-component and the spatial-1 component are 
antisymmetric. Again, together with the parity transformations of the supercurrents, this 
gives the correct parity behavior of the gauge fields, see discussion in Appendix In 
particular, the requirement of a well-defined parity leads to the condition j = 0. 

Inserting Eqs. ( [4.23| ) and the corresponding field strengths into the free energy ( 3.13| ) 
yields 



3 

with Qf, given in Eq. ( 3.14| ) and 



if - /i?) p{B, B) - 2fiB[l^i 7]+{B, B) + j7?_(S, S)]} + , (4.26) 




r]±{B, B) = 2 sinh — — — ^ . ^_ . (4.27b) 
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Again, the asymptotic values of these functions are given in Table |2| in Appendix | 
Minimization of with respect to j yields 

and the minimum of the free energy becomes 



(4.28) 



17 71 



p{B,B) 



(4.29) 



We see that for vanishing magnetic fields the free energy depends on the isospin chemical 
potential, giving rise to a nonzero isospin density. This is expected from the quark content 
of the charged pion condensate and was also observed within the Sakai-Sugimoto model in 



Ref. |6|. 

Taking into account the Meissner effect, which is enforced by the infrared divergence 
in as discussed at the beginning of this subsection, we have to set B = B = 0, leading 
to the simple result for the free energy 



(4.30) 



From Eq. ( ^^Sj) we conclude that for B = B = we have i = 0. And, from the definitions 
(4.17), we see that in the absence of a magnetic field also the normal currents vanish, 
J = J = In the following we shall discuss the results of the pion phase only in the 
presence of the Meissner effect. 

4.3 Meson supercurrents and number densities 

We have seen that all currents in the charged pion phase vanish (except for the supercur- 
rents in the transverse 1- and 2-directions which cancel the applied magnetic field). The 
supercurrents and normal currents in the sigma phase, given in Eqs. (P^ ) and (|1|), 
respectively, are shown in the left panel of Fig. § as a function of the magnetic field. We 
have used the electromagnetic field .Bem, defined in Eq. (3^), with the electric charges of 
up and down quarks. We see that the supercurrents behave linear in i?cm for small i?em 
and approach an asymptotic value for a large magnetic field. These limit cases assume very 
simple forms in terms of the electric quark charges qi and the quark chemical potentials 
A*i,2 = /^B i /^/- The quark supercurrents ji^2 = i i J in the sigma phase then are 



1 

— - < 
2 



3 

Mi sgn Qi 



for small -Bem 
for large Bem 



(4.31) 



The limit of a large magnetic field is strictly speaking not consistent with our approxima- 
tions. Firstly, we have expanded the DBI action for small gauge fields. Secondly, we have 
treated the flavor branes as probe branes which becomes questionable for large magnetic 
fields since one would have to consider the backreaction on the background geometry. As 



- 22 - 



eBe 



4 

eBp, 



Figure 2: Left panel: Meson supercurrents f (red solid) and 2" (blue solid) and normal currents 
J" (red dashed) and J" (blue dashed) as a function of the dimensionless magnetic field eScm in 
the sigma phase. The units are Myj^ja for J' , f and AtM^j^/a for J'^ , J" . Right panel: baryon 
and isospin number densities as a function of the dimensionless magnetic field in the sigma phase. 



The analytical expressions for the functions are given in Eqs. (4.15), (4.18), and (4.33). We have 
fixed [iB = 2/^7 = Mj^^a. In the (charged) pion phase, all currents as well as the baryon density 



vanish due to the Meissner effect; the isospin density is given by the simple expression (4.36b) 



mentioned in Ref. |^] in the same context, the case of a large magnetic field within the 
present approach can only be meaningful if one thinks of the action ( |2.11 ) as a "bottom- 
up" model for QCD, which is not derived from an AdS/CFT correspondence. Indeed, with 
appropriate functions k{z), h{z), the bottom-up model from Ref. can be recovered 
from Eq. (|3.§| ). Thus, in the following we shall use our analytical functions to discuss the 
whole range of magnetic fields with this qualification in mind. 

We can compute the baryon and isospin densities from the free energies computed in 
the previous sections via 

(4.32) 



riBJ 



We obtain for the sigma phase 



AB+Bf , AB-Bf 

{t^B + + {^^B - f^i)- 

p+ p~ 

AB+Bf , AB-Bf 

[fJ-B + IJ-I) [fJ-B - Pi) 



p+ 



p- 



(4.33a) 
(4.33b) 



which agrees with the 0-component of the four-current defined in Eq. ( 4.1(]| ), 



T3 



(4.34) 



provided we use 2/^^,7 (as opposed to ij,b,i) for the boundary values of Aq, Aq, see Table |. 
Remarkably, with these simple rescalings we can reconcile the thermodynamic interpreta- 
tion of as a free energy with the standard definition of the currents from gauge/gravity 
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duality. The reason for this unconventional normalization appears to be the anomalous 
nature of baryon and isospin number in the present model'^. 

As for the 3-components of the gauge fields ( 4.20| ), also the next-to- leading terms in 



the asymptotic expansions of the 0-components only contain the Yang-Mills contribution 
J% ofEq. ( |4l7i| ). 



The densities are plotted in the right panel of Fig. ^. As expected, both densities vanish 
in the case of a vanishing magnetic field. Switching on a magnetic field induces currents 
as well as nonzero densities. Again, it is convenient to express the number densities in 
terms of the quark flavor components, ni^2 = i nj, rather than in baryon and isospin 
components. We obtain for small and large magnetic fields 



^I^h^i^SL forsmahSem 

3 . (4.35) 

fJ-i\qi\Bem for large B 

em 



For the pion phase we find from the free energy ( [4.30| ) 

nl = 0, (4.36a) 

n? = i^^,. (4.36b) 
vr 

From the baryon and isospin densities we can immediately deduce the electric charge 
density nq = qiui + q2n2- The electric charge of the system is relevant for example 
in the astrophysical context because in a neutron star the overall electric charge has to 
vanish. Here we simply observe which electric charge is carried by our system for given 
chemical potentials. For more realistic applications one would have to require charge 
neutrality and possibly counterbalance the charge of the chiral condensate for instance 
by the presence of electrons or protons. For the a phase we find Uq = for vanishing 
magnetic fields, as expected. Switching on a magnetic field induces electric charges in the 
system. For infinitesimally small Sgm a straight line = — 9/i//7 appears in the /i^- 
/i/ plane dividing the plane into a region with infinitesimally positive (above/right of the 
line) and negative (below/left of the line) charge. With increasing magnetic field, giving 
rise to larger charges, the slope of the line slightly decreases and approaches the value 
Hb = —5iJ,i/3 asymptotically for large -Bem- For the pion phase we have Uq = n], which is 
positive (negative) for positive (negative) isospin chemical potentials and independent of 
the baryon chemical potential. 



We may finally recover the scenario considered in Ref. [50| as a limit of our more 



general results. In that paper, a vanishing isospin chemical potential, a vanishing baryon 



This issue is in fact related to the remark below Eq. (4.19): using a modified action as in Ref. |49|] 
which drops certain surface terms would allow us to use the more natural boundary values /is,/- We have 
checked that, apart from the difi'erence in the axial current, the use of the modified action does not lead 
to qualitative changes in our following results. More precisely, changing the action by hand as in Ref. 
amounts to multiplying the Chern-Simons action by 3/2, and rescaling the chemical potential by a factor 



1/2. For the phase diagram presented in Sec. 4.4, only the latter modification is essential and leads to 
corresponding quantitative changes. 
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component of the magnetic field, and an isospin magnetic field constant in the holographic 



coordinate z (as in Eq. (4.22)) was considered. For a comparison it is thus instructive 
to compute the energy density e of the sigma phase. We write the free energy ( 4.13| ) as 
O = e — iiB^B — Pi^i with ub, nj given in Eqs. ( 4.33| ). Then we can express the energy 
density in terms of the number densities, 



32kM2j^ 



a\2 



{B + BY 

For small and large magnetic fields we obtain 



+ (n|-n?) 



{B - BY 



(4.37) 



3iV. 



37r2 

Nr 



{B + BY {B-BY 
(n^ + n^Y {n% - 



\B + B\ 



+ 



\B-B\ 



for small B, B 
for large B, B 



(4.38) 



where we have reinstated the dimensionful magnetic fields according to Eq. (4.10). For 
large magnetic fields we thus obtain, up to a numerical prefactor, an equation of state as 
for a free fermion gas in a magnetic field: setting = i3 = we have Eo- = 6tt'^ {n'^Y / ^c) 
while for a free gas eo = 'n'^rY^/ {BNc) |5^]. In Ref. [^] even the prefactor is exactly that 
of the free gas. We have checked that this discrepancy comes from the surface term in 
the Chern-Simons action: had we dropped the contribution of the surface term (|3.11cj ) 
for the free energy, we would have reproduced Eqs. (34) and (35) of Ref. exactly 
(upon setting n" = B = Q in our result). We can also explain this discrepancy in another 
way. We have computed the free energy, taking into account the Yang-Mills and Chern- 
Simons contributions, and then computed the baryon density (and isospin density) by the 
thermodynamic relation ( [4.32]) . Given this definition of we can write 



''B = §^{jB + 3B) 



967r2 



dz eAfAfPQTr[j^MAf-^PQ] 



(4.39) 



where M, N,P,Q = z, 1, 2, 3, where the first equality can be read off from Eq. ( [4.13| ) upon 
reinstating the dimensionful magnetic field, and the second equality follows from the field 
strengths in Eqs. ( |B.S| ). The resulting expression on the right-hand side differs by a factor 
1/3 from the one in Ref. |5^]; the latter is normalized such that the baryon number is 
equal to the instanton number for a static instanton configuration fl^ (in the absence of 
a magnetic field). Consequently, had we required this normalization we would have had 
to adjust the baryon chemical potential to be S^b with hb the baryon chemical potential 
used above. In this case we would have reproduced the small B limit for the energy density 
of Ref. |5C] exactly. However, in the large B limit we would then have obtained = 2eo/3. 



The reason is that, due to the surface term, the functional dependence of the energy density 



on the magnetic field differs from Ref. |5C], deviating by different factors in the two limits 
of small and large magnetic fields. 
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4.4 Phase diagram and critical magnetic field 

In this section we determine which of the two phases is favored for given values of the chem- 
ical potentials and the magnetic field. To this end, notice that the holographic description 
of our system is our "microscopic" theory; therefore, we have identified the boundary val- 
ues of the bulk field strengths as Bern, not as TCcm- We are interested, however, in a free 
energy comparison at a fixed external magnetic field TCcm- Consequently, we have to apply 
a Legendre transformation to construct the Gibbs free energy G. (See for instance Ref. |15| 
for an analogous construction.) In the case of the charged pion condensate, where Bem = 0, 
the Gibbs free energy is identical to the above computed free energy. 



It is convenient to introduce dimensionless free energies uJa,7r via 



Kk 



KLOn 



(4.40) 



(4.41) 



As we shall see below, k is the only parameter of the model on which the structure of the 



phase diagram depends, see Eq. (4.52). The other constants Mkk and a only set the energy 
scale. To make this k dependence explicit, we have pulled the dimensionless constant k, 
out of uJa,TT- The dimensionless free energies are, using Eqs. (|l2ll) and (|I30| ), 

TT 



{b + bY 



+ (/i 



B - fJ-l) 



{B-Bf 



P+{B,B) P-{B,B) 

where we have introduced the dimensionless chemical potentials 

a/is ~ _ Ci^j 



(4.42a) 
(4.42b) 



fJ-B 



PI 



(4.43) 



Mkk Mkk 

To obtain the Gibbs free energy in the sigma phase we add the contribution B'^^/2 and 



Legendre transform the free energy with respect to the change of variable Be 
Consequently, 



Gcr 



1. 



7^pm ^o" Bq^TC 



For Si given externsil field. TYem 

one determines Bem 
dG^ 



em ' '•em • 

from the stationarity condition 







dB, 



B, 



Ma -He 



where we defined the magnetization in the sigma phase 



Ma 



dBe 



a 



KMa 



Here, the dimensionless magnetization is given by 



Ma 



QiPi 



B + B ( B + Bdp 



P+ 



1 - 



2p+ 




+ <12P2 



B-B ( B-Bdp- 



2p- dI3 



He 



(4.44) 



(4.45) 



(4.46) 



(4.47) 
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where we have used Eq. (|3.2D for the electromagnetic field, and expressed the derivatives 
with respect to B through derivatives with respect to B. Before coming back to the Gibbs 
free energy let us discuss the magnetization and the resulting magnetic susceptibility Xa- 
To obtain M^- as a function of the external magnetic field we first solve Eq. ( 4.45| ), which, 
in dimensionless quantities reads 

i^em = 5cm-«M,, (4.48) 

numerically for i?cm- (Here, the dimensionless field is defined analogously to the field 
i?em, see Eq. ( [4.101) .) Then, we insert the solution back into Eq. ( 4.47 ). The result depends 



on K, for which we have to choose a numerical value. In order to get some numerical 
estimates from our following results we also need to assign values to the other parameters 
of the model. Following 18, Q we shall use 

K = 0.00745 , Mkk = 949 MeV , (4.49) 

which has been obtained from fits to the rho meson mass and the pion decay constant. 
From this value for k we obtain, with Nc = 3 and Eq. ( p. 10 ), the 't Hooft coupling A ~ 16.6, 
and then, with Eq. (|3.7l), a « 2.55. 



The full numerical result for the magnetization is shown in Fig. |^. Our result is in 



qualitative agreement with Ref. |49| , where the magnetization was computed for a one-flavor 
system. (Note, however, that in this reference the boundary value of the field strength was 
interpreted as H, not B.) We see that the magnetization behaves linearly for small magnetic 
fields. The slope is the magnetic susceptibility, i.e., ^ XaHcm- Upon expanding Eq. 



( 4.47 ) for small magnetic fields we find 



'''' 9 l-2|E(,2^2 + ^2^2)- i4.50) 

Since we neither expect the susceptibility to diverge nor to change sign, this result can 
only be trusted for sufficiently small chemical potentials, roughly speaking jlf <C l/(«:e^). 
Given the numerical value (/te^)"^/^ ~ 38 and given that one unit of the quark chemical 
potential /Ij = 1 corresponds to fii ~ 400 MeV, this is not a severe restriction for realistic 
values of /ij. However, this result shows that in principle one has to be careful with large 
chemical potentials in the present approximation where we not only have expanded the 
DBI action for small gauge fields but also neglected the backreaction of the branes to the 
background geometry. 

For large magnetic fields the magnetization saturates. From Eq. ( |4.47|) we find the 
constant value 

Ito M,= «5i-«M, (4.51) 

where we have used that B > B for a two-flavor system with up and down quarks. 

We now return to the free energy comparison. For the sigma phase we insert Eq. ( 4.45| ) 



into the Gibbs free energy (|1|). Then we obtain the following difference in Gibbs free 
energies 

-^''-^^ - -\Hl^ + \>^Ml + n{u^-u:^). (4.52) 



<k/«' - MiJa^ 2 2 
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Figure 3: Dimensionless magnetization Ma- for the sigma phase as a function of the dimensionless 
magnetic field iJem for two different values of the isospin chemical potential and a baryon chemical 



potential fls = 2. The dashed lines are the susceptibilities x<t from Eq. (4.5C), which approximate 



the magnetization for small magnetic fields, Ma ~ Xo-^cm, and the asymptotic values from Eq. 



(MiD. 



If AG > (AG < 0) the vr (a) phase is preferred. It is instructive to relate the comparison 
between the sigma and pion phases to a usual superconducting material, say a metal, where 
we compare the superconducting phase (corresponding to the pion phase) and the normal- 
conducting phase (corresponding to the sigma phase). With the help of this analogy we 
can understand the various terms in AG. The term quadratic in -ffcm is negative, i.e., it 
works in favor of the normal-conducting phase. This term is the free energy cost which the 
superconducting phase has to pay for creating a counter magnetic field in order to expel 
the external magnetic field. In a usual superconductor, this term is thus responsible for 
a critical magnetic field beyond which the Cooper pair condensate breaks down. There 
is an additional term, working against the normal-conducting phase, proportional to M^. 
This term is absent in most usual superconductors which, to a good approximation, have 
no magnetic properties in their normal-conducting phase. We thus expect a competition 
between the two terms, i.e., between the costs that the sigma and charged pion phase 
have to pay for the magnetization and the Meissner effect, respectively. This competition, 
together with the difference iVa — uj-k, will determine the resulting phase diagram. 

For small magnetic fields, -f^cm ^ ^-nd dimensionless chemical potentials of order 
one, i1b,i < 1, we can discuss the phase transition between the sigma and the charged 
pion phase analytically. In this case, the term k^M"^ is of the order of K^e^ and thus 
negligible. The free energy of the sigma phase noja is proportional to ne^ and thus also 
small compared to the remaining terms. We are left with the simple result 

AG , 8k .2 



At the phase transition AG = we thus have 



HL + — • (4.53) 



ixi = ±aP^ ^ ±5.13 /7em . (4.54) 
V K 4 
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Figure 4: Critical magnetic field for the phase transition from the pion to the sigma phase as a 
function of the isospin chemical potential for baryon chemical potentials — 4 (red curve) and 
/is = 40 (blue curve). The dashed line is the approximation from Eq. ( 4.54 ) and almost coincides 
with the curve for /is = 4. Our model does not include a finite pion mass. It can be expected that 
the effect of the pion mass shifts the curves such that they start at fij — dzm,^^ (corresponding to 
jii ~ ±0.375 with the choice ( 4.49D ) instead of at fij — 0. 



This relation can be read as an equation for the critical magnetic field for a given isospin 
chemical potential, or as an equation for the critical chemical potential for a given magnetic 
field. We see that, in this approximation, the phase transition is independent of the baryon 
chemical potential. In Fig. ^ we plot the critical magnetic field as a function of /x/ for two 
different values of /i^- 

Let us now discuss the resulting phase diagram in the plane. Firstly, we consider 

the case of a vanishing magnetic field. Hem = 0. From Eq. (4.53) we see that, in this case, 
the pion phase is favored in the entire ns-fJ-i plane except for the fiB axis. To understand 
this result we recall several features of our model. We treat the fermions as massless (in 
most applications of the Sakai-Sugimoto model, this approximation is used; for a discussion 
about how to incorporate finite mass effects into the model see Ref. |85]). Therefore, 
a charged pion condensate appears for arbitrarily small isospin chemical potentials and 
not only beyond a finite threshold given by the pion mass. Moreover, since we consider 
the confined phase, we cannot account for phases where there is a vanishing (uu) and a 
nonvanishing (dd) condensate or vice versa. In other words, we cannot connect the up- 
flavor branes and leave the down-flavor branes disconnected, as done in Ref. |6^, where 
the deconfined (but chirally broken) phase was considered. And finally, in our approach we 
do not see a phase transition to the chirally restored phase. Since we are in the confined 
phase, where the subspace of the compactified extra dimension X4 and the holographic 
coordinate z is cigar-shaped, the D8 and D8 branes must connect, i.e., chiral symmetry 
must be broken for all values of the chemical potentials. The chiral symmetry can only be 
restored above the deconfinement phase transition at Tc = MKK/(27r), i.e., using ( 4.4S| ), 
Tc ~ 150 MeV. Taking into account these restrictions, our phase diagram at vanishing 
magnetic field is in accordance for instance with Refs. |65, 79, 
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Figure 5: Phase diagram for tiie a and tt phases in the plane of baryon and isospin chemical 
potentials. We have chosen two different values of the dimensionless magnetic field, i/cm = 0.5 
(solid lines, dark shaded sigma phase) and Hem — 5 (dashed lines, light shaded sigma phase). The 
sigma phase contains meson supercurrents while the pion phase contains an isotropic tt^ condensate 
and exhibits the Meissner effect. All lines indicate first order phase transitions. The units of this 
plot, upon fitting the parameters of the model according to ( 4.49 ), are jiBj ~ AiB,//(370MeV) and 
Hem ~ 'Wcm/(l-8 • 10"'^^ G). Hence, due to the huge scales, even compared to magnetar scales, this 
phase diagram is rather of academic interest; for more realistic chemical potentials and magnetic 



fields, the simple approximation (4.54) is sufficient 



Next, we discuss the case of a nonzero magnetic field. The phase diagram for two 
different magnetic fields is shown in Fig. |5|. From Eq. ( 4.54 ) we see that a region for 
the sigma phase opens up, with straight phase transition lines independent of fiB- These 
lines start to bend for larger magnetic fields. We may use the numerical values of the 
parameters of the model given in Eq. ( 4.49| ) and below that equation for some (very) rough 
quantitative predictions from this phase diagram. First we notice that a dimensionless 
field -ffem = 1 corresponds to^ Hem — 2 • 10^^ G, about 4-5 orders of magnitude larger 
than the surface field of magnetars, and most likely even larger than the magnetic field in 
the interior of the star. For the chemical potentials we find that jlsj = 1 corresponds to 
Hbj — 400 MeV. As a comparison, a typical baryon chemical potential for neutron stars is 
at most fiB ^ 1500 MeV, corresponding to flB — 4. Now, as a rough estimate, let us assume 
an isospin chemical potential of 1/10 times the baryon chemical potential in a neutron star, 
i.e., flj ~ 0.4. Then, the phase transition from the charged pion phase to the sigma phase 
occurs at a very large magnetic field of approximately Tiem — 1-6 • 10^^ G. In other words, 
the charged pion condensate at finite isospin chemical potentials appears to be, in terms 
of realistic values for the magnetic field, very robust. 

The superconducting properties of a charged pion condensate have been studied in 
conventional chiral models [37, p^, 8£, 9C]. There, for an isotropic charged pion condensate, 



the scale of the critical magnetic field is set by m^, which is of the order of 10^^ G, whereas 
an anisotropic charged pion condensate has been argued to give a critical magnetic field of 



^In the natural Heavyside-Lorentz system of units of particle physics, a magnetic field strength of 1 eV^ 
corresponds to 51.189. . . G in the Gaussian system. 



-30- 



the order of 10^^ G. In our model, in the vr phase we have an isotropic vr condensate and 
a vanishing 111^^ (see in this context Ref. but nevertheless we have obtained a critical 

field of comparable magnitude. It should be noted that in conventional chiral models the 
charged pion condensate has been found to behave as a type-II superconductor [37, 88, |89|, 
which means that there is another, smaller critical field strength, above which the magnetic 
field can penetrate in the form of magnetic flux tubes. By considering only homogeneous 
fields, we have of course not taken this possibility into account. 



5. Summary and outlook 

We have employed the Sakai-Sugimoto model to study chiral symmetry breaking in a 
two-fiavor system with a magnetic field and baryon and isospin chemical potentials. The 
model consists of a D4/D8-D8 system in which a bulk gauge symmetry on the D8 and 
D8 branes corresponds to a global (flavor) symmetry at the boundary which is interpreted 
as chiral symmetry. Here, left- and right-handed fermions are separated by a fifth extra 
dimension. The electromagnetic subgroup of this flavor symmetry group has been used to 
incorporate an external magnetic fleld. Our results are independent of temperature, valid 
for temperatures below the critical temperature for chiral symmetry breaking. We have 
discussed how the model can account for different Goldstone boson condensates, two of 
which we have described within our ansatz of abelian gauge fleld components. Starting 
from the D-brane action, consisting of Yang-Mills plus Chern-Simons contributions, we 
have solved the equations of motion analytically and computed the Gibbs free energies of 
these two phases. 

The first phase, briefly termed sigma phase, is a linear combination of the usual chiral 
sigma condensate and a neutral pion condensate. We have found that for nonzero mag- 
netic fields this phase exhibits nonzero meson supercurrents. In addition, normal currents 
are generated. For both types of currents, there exist counter-propagating currents such 
that the net vector current of the system is zero. This is reminiscent of unconventional 
(anisotropic) superfluids and superconductors in condensed matter systems or deconfined 
dense quark matter. We have found that the baryon and isospin densities in the system, 
which obviously vanish without a magnetic field, become nonzero once a magnetic field is 
switched on. As a consequence, also an electric charge appears in the sigma phase. 

In the second phase, briefly termed pion phase, charged pions form a condensate. This 
phase reacts very differently to a magnetic fleld. It acts as an electromagnetic supercon- 
ductor, and thus expels the magnetic fleld due to the Meissner effect. We have seen that 
the assumption of a nonzero magnetic fleld would lead to infrared divergences in the en- 
ergy density which cannot be removed by holographic renormalization. Therefore, we have 
introduced supercurrents which induce a magnetic fleld opposite, but equal in magnitude, 
to the externally applied fleld. Then the total magnetic fleld in the vr phase vanishes, and 
a consistent treatment without divergences is possible. In contrast to the a phase, the vr 
phase, due to the Meissner effect, appears unaltered under the influence of a magnetic fleld. 
In particular, the baryon number and all currents (except for the supercurrents cancelling 
the external magnetic fleld) vanish. 
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Besides the calculation of the supercurrents and the observation of the Meissner effect 
in the charged vr condensate, the main result of our work is the free energy comparison 
between the two phases and the resulting phase diagram in the hb-I^i plane. For a vanishing 
magnetic field, a nonzero isospin chemical potential leads to the rotation of the sigma 
condensate into a charged pion condensate. This is expected from studies using the same 



and other models 77, 79, 91, 92 1. In the present study, which does not include quark 
masses, this means that in the absence of a magnetic field the pion condensate is favored 
over the sigma condensate in the entire hb-I^i plane. For a nonzero magnetic field the 
rotation is partially undone, i.e., for a given external magnetic field, there is a region for 
sufficiently small /i/ where the a phase is favored over the vr phase. This is not unlike the 
transition in a metal from its superconducting to its normal conducting state. We have 
found that for small magnetic fields, the critical magnetic field for this phase transition 
is linear in m and independent oi fiB, He oc As a quantitative estimate from our 

result we have discussed that for magnetic fields on compact star scales, the charged pion 
phase at nonzero fij is very robust. More precisely, magnetic fields of the order of 10^^ G 
(well beyond surface magnetic fields of magnetars, although conceivable for their interiors) 
are needed to induce a phase transition from the vr to the a phase for isospin chemical 
potentials of the order of 150 MeV. 

There are several interesting extensions to our work. One could consider the case of 
not maximally separated D8 and D8 branes. This would allow for the possibility of the 
chirally restored phase even at small temperatures and thus might answer the question of 
how a magnetic field affects this symmetry restoration. This extension requires to solve 
an additional equation of motion for the geometry of the D8 and D8 branes, probably 
making an analytical solution impossible. One may also study the question whether the 
solutions found here (anisotropic but homogeneous) are stable against formation of crys- 
talline structures. Moreover, since charged pion condensates have been found to behave 



as type-II superconductors in conventional chiral models [g^, 88, 89|, it would be inter- 
esting to consider inhomogeneous vortex- like configurations of magnetic fields. It is also 
important to check whether the states we have described are stable with respect to other 
meson condensates. We already know that without magnetic field a rho meson condensate 
is expected to form for sufficiently large isospin chemical potentials Another problem, 
not directly related to the present study, but intimately related to the interplay of chirality 
and a magnetic field, is the chiral magnetic effect [|^]. It seems that the Sakai-Sugimoto 
model provides all ingredients to study this effect. 
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A. General form of equations of motion 



Here we present the general form of the variations of the Yang-Mills Lagrangian £ym and 
the Chern-Simons Lagrangian £cs (which can be read off from Eqs. ( |2.12| ) and ( 2.161 ) ) for 
the chirally broken phase. The variation of the Yang-Mills Lagrangian is 



2nVi{2Tia'fTi {tad, - i[ta^ Ap])e-'^^g^'^gP"Tc 



CT/J, 



(A.l) 



where the greek indices run over 0, 1,2,3, u, and where to = 1/2, ta = according to 

the convention ( |2.14 ). Consequently, 



du 
du 
do 



5/2 ^1/2 



V 

5/2 jl/2 



V 



Fm 1 + do 



-ul/2/1/2 
V 



^iO 



Foi + d. 



^ 7,7 



(A.2a) 
(A.2c) 



and 



S'Ukk '^^ym 



^5/2 jl/2 



R^i 



{Sacdu + Aleabc) ^ F„'o + {^acd^ + A'^Cabc) ^y^ F[^ , (A. 



3a) 



.5/2/1/ 



R^t 



{Sacdu + A^eabc) -^ui + {^acdo + ^O^abc) ^^^2/1/2 

+ (5.e9,+46,,,)-j^i^,. 



Oi 



(A.3b) 



{^acdo + Aleabc) 



^5/2jl/2 



F^^ + {5acdi + Al 

^abc) 



^5/2^1/2 



FL , (A.3c) 



where the indices i,j, k run over 1,2,3, and where we used 



KK 



3n 



3/2 
KK 



(A.4) 



with K defined in Eq. ( 3.10 ). 

The variations of the Chern-Simons Lagrangian with respect to the 1 and T3 compo- 
nents are 



5C 



cs 



6A, 
6A- 



.Ka 
~2 



(A.5a) 
(A.5b) 
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with a defined in Eq. (3^). Consequently, 



and 



Mo 

5A'i 
6A?. 



iKa{F:,F^k + FuiFjk)e'"' , (A.6a) 
ma (2F^oF:k - KoF-k + '^FjoFuk - KoF.k) e'^'' , (A.6b) 
iKa{FS,Ff^ + FoFjk)e'''' , (A.6c) 



iKaiF^iFjk + Ff^Fui)e'^>' , (A.7a) 
iKa (2F4F,o - i^^^«o + 2Ff,F.,k - F^^Fju) e^^^ , (A.7b) 
ma{FS>F,k + F%F,o)e'^^ . (A.7c) 



As mentioned in the main part of the paper, we consider maximahy separated branes 
L = tt/Mkki for which the embedding of the D8 branes is trivial, duX4 = 0, and thus v = 1 
(see Eq. ( ^.lOP ). This simplifies the above expressions and also ensures that there is no 
equation of motion for 2:4 (n). The expressions (|A.2|) , ( [A.3|) , ([A.6D , (|A.7|) are used in Sec. 



3.1 to derive the field equations for our specific ansatz. 



B. Solving the equations of motion for constant magnetic fields 

In this appendix we solve the equations of motion for a constant magnetic field for general 
boundary conditions. The resulting general expressions are instructive to see the structure 
and symmetries of the solution. By inserting the boundary conditions from Table |^ into 
the general expressions we obtain the solution for the sigma phase, see Eqs. ( 4.11| ) (the 



pion phase requires a nonconstant magnetic field and is discussed in Appendix The 
general boundary conditions used here are denoted by 

io(±oo) = ^ , ^o(±oo) = 2/if , (B.la) 
is (±00) = 2j^'^ , ^3 (±00) = 2/'^ , (B.lb) 

where the upper (lower) sign corresponds to R (L). It is convenient to express the boundary 
values in terms of their sums and differences, 



V _ f^B,i + f^B,i A _ 
f^B,l = 2 ' ^-S'^ ^ 2 ' (B.2a) 

_|_ ■^L ^R I ^R 

Here, V and A stand for the vector and axial parts of the chemical potentials. 
The general solution for (4.8) with the magnetic field ( [4. 91 ) is 



F^ = ciC+' + C2C+ , F^ = diC' + d2C- , (B.3a) 
F+ = -ciC+' + C2C+ , F^ = -diC^ + d2C- , (B.3b) 
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with constants ci, C2, di, d2 and with 



C±{z) = e^^^^'^ ^ . (B.4) 
Consequently, from Eqs. (|4.7|) we obtain 



kF^o = ciC+^ + C2C+ + diC^ + f^2C- , (B.5a) 

kF,o = ciC+' + C2C+ - diC' - d2C^ , (B.5b) 

feF.s = -ciC+^ + C2C+ - diC' + d2C- , (B.5c) 

A:F.3 = -ciC+' + C2C+ + diC' - c^2C- • (B.5d) 

Here and in the remainder of this and the fohowing appendices we often omit the argument 
z in the various functions for the sake of brevity. For the integration of the field strengths 
we use 

r,,^^.m. ,B.6) 

J k{z) B±B J k{z) B±B 
This yields the gauge fields 

Aq = — X — \- h h ao , (B.Ya) 

B+B B+B B-B B-B 

A giC+^ ^ C2C+ ^ ^iC^ d2C- ^ „ . 

^0 = — - — K h h ao , (B.7b) 

B+B B+B B-B B-B 



^3 = + h h h as , (B.7c 

5 + B 5 + 5 B-B B-B 



^3 = ^^-^ h ^ h as , (B.7d) 

B + B B + B B-B B-B 



with integration constants ao, ao, 03, a^. We determine the eight constants from the eight 
boundary conditions (|B.1|) . This yields the gauge fields 

io = 2ij^ + fii{S+ + S^) + fif{S+ - S„) 

+i(c+ + c_ - r+) + j(c+ - c_ - r_) , (B.sa) 

Ao = 2fiY + fif{S+ + S^) + fii{S+ - S^) 

+ j{C+ + C_ - r+) + j(C+ - C_ - T_) , (B.8b) 
A3 = 2J + j(S+ + 5_)+j(S+-5_) 

+ + c_ - r+) + - c_ - r_) , (b.sc) 
A3 = 2 J + + 5_) + - 

+ 4{C+ + C_ - T+) + - C_ - T_) , (B.8d) 
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and {k{z) times) the field strengths 



kF,o = ^^i[B{C+ + C_) + B{C+ - C_)] + ^^f[B{C+ - C_) + B{C+ + 

+ j[B{S+ + 5_) + B{S+ - + iB{S+ - 5_) + B{S+ + , (B.9a) 
kF,o = i^i[B{C+ - C_) + B{C+ + C_)] + ^^f[B{C+ + C_) + i3(C+ - C^)] 

+ j[B{S+ - S-) + B{S+ + + j[B{S+ + 5_) + 5(5+ - , (B.9b) 
feF,3 = 3[B{C+ + C_) + B{C+ - + - C_) + + 

+ + 5_) + B{S+ - S-)] + ^^[5(5+ - S_) + 5(5+ + , (B.9c) 

A:F,3 = 3[B{C+ - C^) + B{C+ + + j[5(C+ + C_) + 5(C+ - 

+ yLi[B{S+ - S-) + B{S+ + S-)] + ^^[5(5+ + S-) + 5(5+ - 5_)] , (B.9d) 



where the functions C±{z), S±{z), and T+ are defined in Eqs. ( 4.12| ). As it should be 



the gauge fields ( |B.g| ) transform as a vector under a parity transformation. This can be 
seen as follows. A parity transformation is given by {xi,X2-,x^, z) — > (— xi, — X2, — 2:3, — z). 
In particular, the transformation z —z implies a chirality transformation i? — > L since 
the two halves of the D8/D8 branes, namely z > and 2: < 0, correspond to right- and 
left-handed fermions. Consequently, a parity transformation acts as C+(z) +(^±(-2), 
S±{,z) — > — 5+(2;) (since the magnetic fields 5, B are even under parity). For the super- 
currents we have +j, +j and J,J^ — J, — J. Here we have used that the Goldstone 
bosons are pseudoscalars (for a detailed discussion of the parity of the mesons in the Sakai- 
Sugimoto model see Ref. fl^). As a result we see that the temporal components ( B.8a| ), 



(B.Sb) have even parity, while the spatial components (B.8c), (B.Sd) have odd parity. This 
statement is true for arbitrary values of the currents j, j, J, J. We shall see below that in 
the case of a charged pion condensate the requirement of a well-defined parity results in 
conditions for the supercurrents, see discussion below Eq. ( p.lC| ). 
In order to compute the free energy we note that 

Cl-Sl = . 1 . (B.IO) 

^ ^ sinh2[7r(5±5)/2] 

Therefore, the following combination of field strengths, needed for the free energy, becomes 
independent of z, 

{-Fl - + Fl + Fl) = 2[(i + jf - i^i + t^f n + 



sinh^[7r(5 + 5)/2] 

+ 2[(i - jf - (^^ - f^ff]--^^j^^^, ■ (B.ll) 

smh^[7r(5 - 5)/2] 

Next we use the fact that S± and C± are antisymmetric and symmetric in z, respectively, 
as well as 

C±(oo) =coth^^^^^-^, 5±(oo) = l, (B.12) 
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small B, B 


large \B\, \B\ 






\B\ > \B\ 


\B\ < \B\ 


p± 


6 ■k{B±BY 
vr D 


2\B±B\ 


p 


12 5^2 + ^2 

\ TT 

vr 15 


A\B\ 


m\ + \B\) 


v+ 


ttBB 
3 


2BsgnB 


B sgn B + B sgn B 


V- 


25 


2B 


{\B\ + \B\)sg-aB 



Table 2: Behavior of the functions p, p±, rj±, defined in Eqs. ( 4.14 ), ( 4.27 ) for small and large 
magnetic fields B, B. We have kept relative magnitude and sign of baryon and isospin components 
arbitrary. They can then later be inserted according to the electric charges of the quarks. We show 
the behavior for small magnetic fields up to second order and the behavior for large magnetic fields 
in leading linear order. 



to obtain 



A^kF^Q + AokF,o - A3kF,3 - A^kF, 



23 



8//^(ji? + jB) + 8f,^ {jB + jB) - SJifi^B + ^ifB) - SJ{fij,B + ^ifB) 



+ 4[(/.^ + ^iff- (i + jf] {B + B) coth 



tt{B + B) 



(B.13) 



Inserting Eqs. ( [B.ll| ) and ( |B.13| ) into Eq. ( 3.13| ) yields the free energy 



(Pb - Pi 



A\21 



~4^lUJB+JB)-4^^Y{jB + JB) 

+ 4Ji^^iB + nfB) + AJ{^iB + /.^^)} , 



(B.14) 



with p± defined in Eq. (4.14). For the behavior of p± for small and large magnetic fields 
see Table |2[ We see that if we allowed for nonzero axial chemical potentials ^f, 
the free energy would be unbounded from below in the directions of the sum of left- and 
right-handed supercurrents J and J. However, in the physical case of the a phase where 
P'b ~ pf ~ ^ ^'■^^ energy remains bounded and independent of J and J. The latter is 
a manifestation of a residual gauge symmetry ( "residual" since we have already employed 
the gauge = 0), i.e., we can choose a gauge where J = J = 0. This is in accordance 
with the discussion in Ref. [17|, see in particular Eq. (5.23) in this reference. 
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Minimization of with respect to the currents j, j yields 



fi^s + f,^ B + B fi^s-f,^ B-B 

2 P+ 2 p_ 

^i^ + ^iY B + B fi^-f^Y B-B 



2 p+ 2 p_ ' 

and the minimum of the free energy becomes (with fi^ = fif = 0) 



(B.15a) 
(B.15b) 



^2 {B + B) 

p+ 



V ,.V\2 



+ (Mb - 



iB-B) 



(B.16) 



C. Solving the equations of motion for nonconstant magnetic fields 



In this appendix we present the general solution to the differential equations (4.8) for the 



case of a nonconstant isospin magnetic field given in Eq. ( 4.22| ). The general expressions 



given below reduce to the results for the charged pion phase upon inserting the specific 
boundary conditions from the second row of Table |l]. The general boundary conditions 
considered here are the same as the ones given in Eqs. ( p. 11 ). 

Then, the solution of (|4.8| ) has the same form as given in Eqs. ( |B.3D and (|B.5| ), with 
C,±{z) replaced by 

= g(Bi-^ arctan z) arctan 2 -y"^ 

To obtain the gauge fields we need 

■iz^Jrlp = <3+ (z) , j = -Q- (z) . (C^2b) 



with P±, Q± given in Eqs. ( |4.25|) . We shall denote Q'^ = Q±{+oo), = Q±{—oo), 
= P±{+oo), P^ = P±{—oo), and use P^ = Pj, = Q^. Hence we can express 
the values of P_, Q- at z = ±oo through the values of P+, at z = ^oo. Then, after 
determining the integration constants from the boundary conditions we can write the gauge 
fields as 

io = 2^1 + f^BiS+ + S-) + j{S+ - S.) 

+ fif{C+ + (7_ - f+) + j{C+ - C_ - f_) , (C.3a) 
^0 = 2fiY + fif{S+ + S^) + j{S+ - S-) 

+ fii{C+ + C„ - f+) r_) , (C.3b) 

is = 2 J + j{S+ + 5_) + fif{S+ - 5_) 

+ jiC+ + f+) + f,iiC+ f_) , (C.3c) 

^3 = 2J + X5+ + 5_)+/i^(5+-5_) 

+ i((7+ + r+) + fif{C+ T_) , (C.3d) 
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and the field strengths as 

kF,Q = /x^(c+ + c_) + j(c+ - c_) + ^lf{s+ + s_) + i(s+ - s_) , (C.4a) 

/cF^o = /x:f(c+ + c_)+i(c+-c_) + /i^(s+ + s_)+j(s+-s_), (C.4b) 

A:F^3 = i(c+ + c_) + /u^(c+ - c_) + + s_) + - s_) , (C.4c) 

/cF,3 = j(c+ + c_) + //;^(c+ - c_) + + s_) + - s_) , (C.4d) 









--P- 











where C±, S±^ and T± are defined in Eqs. (|4.24| ), and where 

c+{z) = ' , c_(z) = ^'^^+ > (C.5a) 

(These additional definitions were not necessary in the case of constant magnetic fields, 
since there the integration of the solution could be expressed in terms of the same functions 
as the solution itself.) 

We now have to check the behavior of the gauge fields ( p.3| ) under a parity transforma- 
tion. For the pion phase we have = /i^ = 0. We have to require Aq +Aq, Aq —Aq, 
A^ — > —^3, ^3 — > +^3 (note the additional "twist" for the isospin components originating 
from the isospin rotation explained in Sec. LI). Since C±{z) — > +C±{z), S±{z) — > —S±{z), 
and J, J — > +j, +j and J, J ^ — J, — J under a parity transformation, we have to require 

J = j = 0. (C.6) 

We shall continue with the general solution but have to keep this condition in mind for the 
final result. 

For the free energy we first note that the following combinations are independent of z, 

^-^- + ^^^- = (p^-p,-kqI-q;)' ^+^- + ^-^- = °- ^^-'^ 

Then, we find 

Next we use that c± and s± are symmetric and antisymmetric in z, respectively, and denote 
= c±(oo) = c±(— oo), = s±{oo) = —s±{—oo). Then, 



AokF^o + AokF^o - A^kF,^ - A^kF^^^ 

= 4(4 + st)iplfif + /i^/i^ - jj - jj) + 4(4 - st)ijfil + Jf^Y - Jf^i - J^if) 
+4(4 + cl )[(m^)2 + {4f - if + /)] . (C.9) 

Inserting this into the free energy ( |3.13 ) yields 

n = [f + / - if^if - i^^ff] P - 2i^^Y4 + pY^^i) ry+ - 2if.lj + ^^Yj) ri- 

+ 2J(/xV +m) + 2J(//V +m)} , (c.io) 
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with p and ry-t = s]]] it si given in Eqs. ( [4.27 ); their behavior for small and large magnetic 



fields can be found in Table ^. As in the case of constant magnetic fields discussed in the 
previous appendix, see Eq. ( B.14D , the free energy is unbounded from below without further 



constraints. This can be seen by computing the matrix of second derivatives d'^Q/{dxmdxn) 
with Xm, Xn = J, J- This matrix has eigenvalues 2kM^^/3 [p ± (p^ + 4r/^)^/^], two of 
which are negative for all magnetic fields. However, we already know from the requirement 
of a well-defined parity of the gauge fields that J = j = 0. Then, with p^ = pj = 0, as 
required for the charged pion condensate, we see that the free energy becomes bounded 
from below. The only remaining supercurrent with respect to which we need to minimize 
the free energy is then j. The sum of left- and right-handed supercurrents, J, remains 
undetermined, which is, as mentioned for the case of the sigma phase below Eq. ( [B.14| ), a 
consequence of the residual gauge freedom. We may thus set J = 0. 
We can now minimize with respect to j, 

J = M^y, (C.ll) 

and insert this back into the free energy, 

^o = -^-^[{f^'^f^^ + {f^f?P + 2v^pW} ■ (C.12) 

D. Equations of motion and free energy in the chirally restored phase 

Within our approximation of treating the fiavor branes as probe branes, the free energies 
discussed in the main part of the paper are negligible for the finite-temperature phase 
transition to the chirally restored phase. It is rather the background geometry which 



is responsible for this phase transition p^, 63]. Therefore, our approach cannot show 



magnetic-field induced corrections beyond the order of Nf /N^ to the critical temperature 
Tc for chiral symmetry breaking. This is different when the D8 and D8 branes are not 



maximally separated in the extra dimension |94, 95 1. 



In this appendix we simply give the equations of motion and the free energy for the 
chirally restored phase without discussing the solutions. We do so for the sake of complete- 
ness but also because these expressions may be useful to compute possible small corrections 
to Tc of the order of Nf/Nc- One might then speculate whether these corrections persist 
for smaller and thus more realistic values of Nc- We leave such a study for the future. 

The derivation of the equations of motion and the free energy of the chirally restored 
phase is analogous to the one for the confined phase given in Sec. |3| and Appendix The 
only difference is the use of the metric ( ^.9b| ) instead of ( 2.9a ) and Eq. (2^) instead of 



( |2.3| ). We use the same coordinate transformation as in the chirally broken phase, i.e., Eq. 
( |3.1| ) with ukk replaced by ut and with z G [0,oo]. This is not really a simplification in 
this case but it helps to compare the result to the one for the chirally broken phase. We 
find for the equations of motion 

d,[h{z)d,b] = d,[k3{z)d,b] = 0, (D.l) 
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and 



dz[ko{z)F^o] 



(27rT)3 
(27rr)3 
(27rr)3 
(27rr)3 



6(z)F,3 + 6(z)i^.3 


(D.2a) 


6(z)F,3 + b{z)F,3 


(D.2b) 


b{z)F,o + b{z)F,o 


(D.2c) 


b{z)F,o + b{z)F,o 


(D.2d) 



In contrast to the confined phase, there are now two different functions appearing for the 
temporal and spatial components, 



ko{z) 



z u 



1/2 



k-i{z) = ZUrl'^iu^ + UTZ^f'^ . 



(D.3) 



The free energy becomes 



deconf odeconf , odeconf 



j^aecont ^ 

2k{2ttT)^ 



+ 



k{2ttT) 



3 /-oo 



dz 



-koiz)iF^o + F^o) + k3iz)iF^3 + 



ko{z){AoF,o + AoF,o) - ksiz^A^F^s + A^F,^) 



z=+oo 

2 = 



where 



deconf 



32K{2TTTf 
9(27rQ')2uf,MKK 



/•oo 

/ dz z (mj. + u^z 
Jo 



j^deconf ^ ><^02^q2^ H dz Z u}'\ul + UtZ^)-''/^ . 
M-KK Jo 



(D.4) 

(D.5a) 
(D.5b) 



Here we have assumed the magnetic field to be constant in z, b{z) = B, b(z) = B, which 
solves Eq. ( |D.1| ). We see that at the critical temperature where 2ttT = Mkk and thus 
ut = ukk the free energy assumes a form very similar to the one in the confined phase. 
The only differences are then the functions kQ{z) and ^3(2) (vs. the single function k(z) in 
the confined phase) and the different integrands in Qg and ^Ib- 
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